Phase-fluctuation induced reduction of the kinetic energy at the superconducting 

transition 
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Recent reflectivity measurements indicated a possible violation of the in-plane optical integral in 
the underdoped high- Tc compound Bi2Sr2CaCu20s+s up to frequencies much higher than expected 
by standard BCS theory. The sum rule violation may be related to a loss of in-plane kinetic energy 
at the superconducting transition. Here, we show that a model based on phase fluctuations of the 
superconducting order parameter introduces a change of the in-plane kinetic energy at T^. The 
change is due to a transition from a phase-incoherent Cooper-pair motion in the pseudogap regime 
above Tc to a phase-coherent motion at Tc. 

PACS numbers: 71.10.Fd, 71.27.+a, 74.25. Jb, 74.72.Hs 



I. INTRODUCTION 

The key idea of the phase-fluctuation scenario in the 
high-Tc superconductors is the notion that the pseudo- 
gap observed in a wide variety of experiments arises from 
phase fluctuations of the superconducting gapLSiMiSiSii. 
In this scenario, below a mean-field temperature scale 
T*^^, a dx2-y2-wave gap amplitude is assumed to de- 
velop. However, the superconducting transition is sup- 
pressed to a considerably lower transition temperature 
Tc by phase fluctuations^'^'. In the intermediate temper- 
ature regime between T*^^ and Tc, phase fluctuations 
of the superconducting order parameter give rise to the 
pseudogap phenomena. Recently, we have shown that 
indeed a two-dimensional BCS-like Hamiltonian with 
a (ia;2_j,2-wave gap and phase fluctuations, which were 
treated by a Monte-Carlo simulation of an XY model, 
yields results which compare very well with scanning tun- 
neling measurements over a wide temperature rangeiSi. 
Thus, they support the phase-fluctuation scenario for the 
pseudogap. 

There is also increasing evidence from a number of 
recent experiments for the relevance of phase fluctua- 
tions such as the Corson et al. measurements® of the 
high-frequency conductivity that track the phase corre- 
lation time. XY vortices are probably responsible for 
the large Nernst efFec1ii2*ii. The evolution of Tc with 
electron irradiation found very recentljJ^ also empha- 
sizes the importance of phase fluctuations. In this pa- 
per, we argue that phase fluctuations should have yet 
another rather unexpected consequence: they induce a 
reduction of the kinetic energy at the superconducting 
transition. This reduction is due to a transition from 
a "disordered", i. e. phase-incoherent Cooper-pair mo- 
tion in the pseudogap regime above Tc to an "ordered" , 
i. e. phase-coherent motion at Tc. Comparison of our re- 
sults, based on the BCS phase-fluctuation model, with 
optical experiment s ^•^1^'* support this idea. 

In ordinary BCS superconductors the optical conduc- 



tivity is suppressed at frequencies within a range of 
about twice the superconducting gap. The correspond- 
ing low-frequency spectral weight Wiow is transfered to 
the zero-frequency delta peak Wn^^, associated with 
the dissipationless transport (and the superfluid weight 
D) in the superconducting state. This is the Glover- 
Ferrell-Tinkham (OFT) sum rule. On the other hand, 
the total frequency integral of the optical conductivity 
is conserved, when decreasing the temperature across 
the superconducting transition, due to the f-sum rulai^, 
i. e. Wi,% = Wp,,. 

However, recent measurements of the in-plane optical 
conductivit}ii2*i^ have indicated a violation of the GFT 
optical sum rule for frequencies up to 2eV in underdoped 
Bi2Sr2CaCu20s+s (Bi2212). By entering the supercon- 
ducting state, not only spectral weight Wiow from the 
microwave and far-infrared, but also from the visible op- 
tical spectrum i. e. high-frequency spectral weight Whigh 
contributes to the superfluid condensate Wd- That is, 
in contrast to ordinary BCS-superconductors, a color 
change" is introduced at the superconducting transition. 
The interpretation of this unusual result may require the 
inclusion of local-fleld effects and other (such as exci- 
tonic) many-body effects. They are known to play a 
crucial role already in weakly-correlated systems (such 
as semiconductors), and introduce a shift of order of 
the Coulomb correlation energy between single-particle 
and two-particle, i. e. optical excitationsi^. Therefore, 
they may partly account for the "high-energy" features 
observed in cr(w). On the other hand, within a tight- 
binding one-band model, the anomalously large energy 
scale, which contributes to the superfluid weight, and the 
corresponding color change can be attributed to a reduc- 
tion of kinetic energyi^, at the superconducting transi- 
tion. This is rather surprising, since one would expect 
that in a conventional (BCS) pairing process, it is the 
potential energy which is reduced at the expense of the 
kinetic energy, with the latter being increased due to 
particle-hole mixing. 

The full optical integral, when integrated over all fre- 
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quencies and energy hands, is proportional to the carrier 
density (n) over the bare mass (m) 



Wd + Whigh = / Re axxii^) duj = -— , 
Jo 2m 

(1) 

and, thus, is conserved. When the optical integral is 
restricted over a finite (low) range of frequencies Q, in 
the HTSC typically of the order of eV, one may consider 
the weight Wiow + Wd as being essentially due to a single 
band around the Fermi energy, i. e. 



Wr 



Re aa;x{(^) duj ^ {ire'^a^ /2h^V)EK, 



(2) 

where a is the single-band conductivity, a the lattice con- 
stant and V the unit cell volume. With this single-band 
assumption, the frequency integral of the optical conduc- 
tivity is proportional to the the inverse mass tensor (^^, 
X being the direction in which the conductivity is mea- 
sured) weighted with the momentum distribution nj^^i^: 



Ek = {2/a^N) 



k 



d^k_ 

dk'- 



-nk, 



(3) 



with N the number of k points. This quantity depends 
upon the bare single-particle band structure e^, being 
proportional to minus the kinetic energy Ek = — -Efcm 
for a (nearest-neighbor) tight-binding (TB) model, while 
for free electrons it is a constant given by the electron 
density divided by the effective mass. 



II. PHASE FLUCTUATION SCENARIO FOR 
KINETIC ENERGY REDUCTION 

In this paper, we propose as a mechanism for a kinetic- 
energy reduction phase- fluctuations. That is, in order 
to have condensation into the superconducting state, 
one needs, in addition to the binding of charge carri- 
ers into Cooper pairs, long-range phase coherence among 
the pairs. Since superconductors with low superconduct- 
ing carrier density (such as the organic and underdoped 
high- Tc superconductors) are characterized by a rela- 
tively small phase stiffness, this implies a significantly 
larger role for phase fluctuations, than in conventional 
superconductorsiiSSiSi. As a consequence, in these ma- 
terials the transition to the superconducting state does 
not display a typical mean-field (BCS) behavior, and 
phase fluctuations, both classical and quantum, may have 
a significant influence on low temperature properties. 
When coherence is lost due to thermal fluctuations of the 
phase at and above the transition temperature Tc, pairing 
remains, together with short-range phase correlations. 
These phase fluctuations can cause the pseudogap phe- 
nomena observed e. g. in tunneling experimentsiSiSSiS 
in the underdoped HTSC. 



We show here that indeed phase fluctuations con- 
tribute to a significant reduction of the in-plane kinetic 
energy upon entering into the superconducting phase be- 
low Tc, with a magnitude comparable to recent experi- 
mental results. The physical reason for this kinetic en- 
ergy lowering is that, due to phase fluctuations and to the 
associated incoherent motion of Cooper pairs (cf. Fig.^l, 
the pseudogap region has a higher kinetic energy than the 
simple BCS mean-field state. When long-ranged phase 
coherence finally develops at Tc , the Cooper-pair motion 
becomes phase coherent and the kinetic energy decreases. 
The onset of the coherent motion can be seen, for exam- 
ple, from the development of coherence peaks in the tun- 
neling spectrum of BiSrCaCuO compounds (see e. gi^ 
and Fig. |2J|. The initial cost of kinetic energy, which 
is needed for pairing, is payed at a mean-field tempera- 
ture T^ ^^ considerably higher than Tc- Therefore, the 
reduction of kinetic energy observed experimentallyi^ii^ 
can be attributed to a transition from a phase-disordered 
pseudogap to a phase-ordered superconducting state. We 
stress that this effect is independent of the particular 
mechanism leading to pair formation, as long as the su- 
perconductor considered is characterized by a small phase 
stiffness^. 

Our starting Hamiltonian is of a simple BCS form given 

by 

H^K-^J2(^^^(^ls)+^ls{^^^))' (4) 

i 5 

with the nearest-neighbor hopping term 

K^-t Y.{cl^cj, + c]^c,^), (5) 

where c]^ creates an electron of spin cr on the z*'* site 
and t denotes an effective nearest-neighbor hopping. The 
(ij) sum is over nearest-neighbor sites of a 2D square lat- 
tice, and, in the pairing term, S connects i to its nearest- 
neighbor sites. The local d-wave gap, 



{Als)^^{4A+Si-4As,)-^e'''-, (6) 

is characterized by the fluctuating phases 

{^Pi + (pi+s)/2 for S in x-direction 



{^Pi + ipi-^-s)/2 + TT for S in y-direction. 



(7) 



and by a spatially constant amplitude A. We neglect 
the relative bond phase fluctuations between S = x and 
y as well as amplitude fluctuations. Thus, we consider 
only center of mass pair phase fluctuations, which are the 
relevant low-energy degrees of freedom, in a situation in 
which the superfluid density is small, like in the under- 
doped cuprates. In our two-dimensional (2D) model Tc 
corresponds to the Kosterlitz-Thouless transition tem- 
perature TfCT, where the phase correlation length ^ di- 
verges. 
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With this Hamiltonian it is straightforward to show, 
that the optical sum rule yieldsi2iS^i2&, 



Re axx{^) du! — —e'^Tr{kx) /2, 



(8) 



in units where h = c — 1, with {k^) being the expec- 
tation value of the nearest neighbor hoppingSi in the x- 
direction, i. e. 



(9) 



Since we are only interested in the temperature region 
T > Tc, one can safely assume that the fluctuations of 
the phase ipi are predominantly determined by a classical 
XY free energy^^^ 



F[ipi] = -jy^cos {cpi 



(10) 



Our physical picture here is that the XY-action arises 
from integrating out the shorter wavelength fermion de- 
grees of freedom, including those responsible for the 
formation of the local pair amplitude and the internal 
dx^-y^ structure of the pair. Thus, the scale of the XY- 
lattice spacing is actually set by the pair coherence length 
^0- In our work, we have chosen A so that £,o ^ ^ I- 
In this case, the phase configurations ifi calculation can 
be carried out on the same L x L {L = 32) lattice that 
is used for the diagonalization of the Hamiltonian. This 
allows the Kosterlitz-Thouless phase correlation length ^ 
to grow over a sufficient range as T approaches Tkt and 
minimizes finite-size effects. Thus, we are always in the 
limit where the phase correlation length ^ is larger than 
the Cooper pair size ^q, when the temperature T is below 
the mean- field critical temperature T*^^. 

In principle, also the coupling energy J can be con- 
sidered as arising from integrating out the high-energy 
degrees of freedom of the underlying microscopic sys- 
tem. Here, we will proceed phenomenologically, neglect- 
ing the temperature dependence of J and simply use 
it to set the Kosterlitz-Thouless transition temperature 
Tkt equal to some fraction of T*^^. Specifically, for 
the present calculations we will set Tkt — jT^^. This 
choice is motivated by the recent scanning-tunneling re- 
sults in Bi2Sr2CuOQ^s, where ~ lOK and the pseu- 
dogap regime extends to about 50K, which we take as 

rpMF 

In a previous paper*" we have presented a detailed nu- 
merical solution of the 2D-BCS like Hamiltonian of Eq.01 
with a d-wave gap and phase ffuctuations. This is a min- 
imal model but, nevertheless, contains the key ideas of 
the cuprate phase fluctuation scenario: that is, a d-wave 
BCS gap amplitude forms below a mean-field temper- 
ature T^^ , but phase fluctuation suppress the actual 
transition to a considerably lower temperature Tc- In 
the intermediate temperature regime between T^^^ and 
Tc, the phase fluctuations of the gap give rise to pseudo- 
gap phenomena. Comparison of these results with recent 
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FIG. 1: Kinetic energy per bond {kx) as a function of temper- 
ature for the non- interacting tight-binding electrons (TB) , the 
BCS solution (BCS), and our phase-fluctuation model (PP) 
for jj, — ((n) = 1). The large vertical arrows indicate the in- 
crease in kinetic energy upon pairing relative to the free tight- 
binding model, and the small arrows indicate the additional 
increase due to phase fluctuations. This additional phase- 
fluctuatton energy rapidly vanishes near Tc = Tkt, which 
causes the signiflcant change in the optical integral upon en- 
tering the superconducting state at Tkt ~ O.lt. Note that 
the thick line follows the actual kinetic energy encountered in 
our model, when going from the pseudogap to the supercon- 
ducting regime. 



scanning tunneling spectra of Bi-based high- Tc cuprates 
supports the idea that the pseudogap behavior observed 
in these experiments can be understood as arising from 
phase fluctuations^. 

In the present calculations, where we assume a BCS 
temperature dependence of the pairing gap A(T), we 
have therefore set A(T = 0) = l.Qt corresponding to 
Tc*^^ ~ 0.42t and selected J so that Tkt = 0.1*22,. The 
condition that ^ > is thus always fulfilled, if we are 
not too close to T^^ . The calculation of the kinetic en- 
ergy for an L X L (T = 32) periodic lattice now proceeds 
as follows22ii22i: a set of phases is generated by a 
Monte Carlo importance sampling procedure, in which 
the probability of a given configuration is proportional 
to exp(-T[(p,]/T) with F given by Eq. With {(^J 

given, the Hamiltonian of Eq. Q is diagonalized and 
the kinetic energy Ekin(T, {ipi}) = {kx){ipi} is extracted. 
Further Monte-Carlo {ipi} configurations are generated 
and an average kinetic energy EkiniT) = (kx), at a given 
temperature, is determined. 

Fig-ffldisplays the kinetic energy {kx) as a function of 
temperature for non-interacting tight-binding electrons, 
for BCS electrons, and for our phase-fluctuation model, 
respectively. We can clearly see that pairing, as expected, 
produces an overall increase of kinetic energy (indicated 
as vertical arrows) with respect to the free-electron case. 
We observe that in the phase-fluctuation model the ki- 
netic energy is further increased (small vertical arrows) 
due to the incoherent motion of the paired electrons. 
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FIG. 2: Kinetic energy contribution from phase fluctuations 
5{kx) = {kx)pp — (kx)BCS- One can clearly see the sharp 
decrease of the kinetic energy near the Kosterlitz-Thouless 
transition at T = Q.lt = Tc. AEk gives a estimate of the 
kinetic condensation energy. 

The kinetic energy is a smoothly decreasing function of 
temperature for T 0. This is expected from the fact 
that, at high temperature, more electrons are transferred 
to higher kinetic energies, and is in agreement with the 
experimental resultsi^*i^. What we are especially inter- 
ested in, is the rather pronounced change (magnified in 
Fig. 121 by using a different scale for the kinetic energy) 
near Tc = Tkt, where the kinetic energy of our phase- 
fluctuation model rather suddenly reduces to the BCS 
value. This sudden deviation from the T > Tc behavior 
is also obtained in experiments, which show a kink in 
the temperature dependence of the low-frequency spec- 
tral weight Wioto + Wd at Tj^. 

This pronounced change of in-plane kinetic energy can 
be better observed in Fig. [3 where we plot the difference 
between the BCS kinetic energy and the kinetic anergy of 
our phase-fluctuation model S{kx) — {kx)pp — {k 

x) BCS- 

As discussed above, this reduction is due to the onset 
of phase coherence of the Cooper pairs below the su- 
perconducting transition temperature Tc = Tkt- This 
is signaled by the appearance of sharp coherence peaks 
in the single-particle spectral function upon developing 
long-range phase coherence^. The corresponding result 
for the density of states 7V(w) displaying these coherence 
peaks is shown in Fig. |21 

Notice that this argument for the reduction of kinetic 
energy at Tc due to a phase ordering transition is quite 
robust. For example, we expect it to be valid (and actu- 
ally to be stronger) in a true three-dimensional system. 
As a matter of fact, it has been arguecii*S£ that even 
small interplane couplings play an important role due to 
the infinite-order nature of the the KT transition. 

In order to get a rough estimate of the kinetic con- 
densation energy, we calculate the reduction in kinetic 
energy near Tc, i. e. 

AE,^-^f {Redll{uj)-RedZA^))du;, (11) 




FIG. 3: Single particle density of states N{u) for different 
temperatures T for a 32 x 32 lattice. Coherence peaks de- 
velop (marked with circles for T = l.STc) as T approaches 
Tc = Tkt in exactly the same temperature regime (Tc < T < 
1.5Tc) where in Fig. the kinetic energy reduction occurs. 
The inset shows the corresponding temperature dependence 
of the phase correlation length ^(T). 

as indicated by the energy change AiJ^ in figure 13 As- 
suming that t ~ 250 meV, we get a condensation energy 
estimate of 1.5meV per Copper site, which is in order 
of magnitude agreement with the experimental results 
(again assuming a one-band TB analysis). 

Up to now, to refrain from further approximations, 
we have set the chemical potential fi equal to zero and 
have only considered nearest-neighbor hopping. We have 
checked to some extent, how robust these results are with 
respect to finite doping ((n) sa 0.9) and the inclusion of 
a next-nearest neighbor-hopping term t' in our Hamil- 
tonian Eq. ^ Notice that in this case Ek is no longer 
proportional to the kinetic energy. For t' < 0.3t, our re- 
sults for the sum rule violation are reduced only by about 
20% - 30%. 



III. SUMMARY 

In conclusion, we have shown that the recently ob- 
served violation of the low-frequency optical sum rule 
in the superconducting state, associated with a reduc- 
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tion of kinetic energy, can be related to the role of phase 
fluctuations. The decrease in kinetic energy is due to 
the sharpening of the quasiparticle peaks close to the su- 
perconducting transition at = Tkt, where the phase 
correlation length ^ diverges. We suggest that this sum 
rule violation should also appear in other superconduc- 
tors with low charge carrier density (phase stiffness) such 
as the organic superconductors. 
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